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Preface 


The  purpose  of  this  thesis  was  to  study  the  effects  of 
turbulent  mean  flow  on  the  frequency  response  of  round  air 
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thesis  advisor,  was  always  helpful  whenever  I  ran  into 
problems.  Harley  Linville  was  always  ready  to  help  when 
equipment  problems  cropped  up.  I  also  wish  to  thank  my 
wife,  Theresa,  for  putting  up  with  me  during  the  final 
months  of  study  and  typing. 
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Abstract 


<■  '^Theoretical  predictions  of  the  small  signal  frequency 
response  of  round  pneumatic  transmission  lines  with 
turbulent  mean  flow  were  compared  with  experimental  results. 
The  frequency  response  curves  were  found  for  lines  varying 
in  length  from  24  to  36  in.  with  inside  diameters  of  0.195, 
0.119,  and  0.041  in.  The  lines  were  tested  at  Reync 
numbers  of  2000,  5000,  and  10000. 

Theoretical  solutions  were  obtained  using  Nich<  > 
equations  as  modified  by  Krishnaiyer  and  Lechner.  Solutlous 
were  also  found  using  several  different  modifications  of  a 
constant  LRC  model  developed  by  Moore.  The  results  were 

mixed:  for  the  0.195  and  0.119  in.  ;  lines  the  prediction  of 

gain  was  good  but  for  the  0.041  in.  lines  the  results  were 

poor.  The  accuracy  and  applicability  of  the  constant  LRC 

model  was  explored  along  with  its  various  modifications. 
The  constant  LRC  model  with  the  AC  resistance  showed 
potential  for  predicting  the  gain  in  fluid  transmission 


lines  with  turbulent  flow. 


The  limitations  and 


applicability  of  the  constant  LRC  models  was  studied 


THE  EFFECTS  OF  MEAN  FLOW  ON  THE  DYNAMIC 
CHARACTERISTICS  OF  FLUID  TRANSMISSION  LINES 


I. Introduction 


Background 

The  propagation  of  small  signals  in  fluid  lines  has 
been  of  interest  to  scientists  and  engineers  for  many  years 
and,  particularly  since  the  advent  of  fluidic  devices  in  the 
1950's.  Many  investigations  have  been  conducted,  both 
analytical  and  experimental,  in  an  effort  to  develop  models 
that  can  accurately  predict  the  frequency  response  of  fluid- 
filled  transmission  lines.  Equations,  develrped  by  Nichols 
(Ref  1)  and  modified  by  Krishnaiyer  and  Lechner  (Ref  2), 
have  proven  useful  in  predicting  the  frequency  response 
under  laminar  flow  conditions. 

Brown,  Margolis,  and  Shah  (Ref  3)  using  two  and  three 
region  viscosity  profiles  predicted  increased  attenuation 
with  Reynolds  number  in  lines  with  fully  developed  turbulent 
flow.  These  models  however  resulted  in  relatively 
complicated  solutions.  Moore  (Ref  4)  used  a  simple 
constant  LRC  model  as  suggested  by  Brown  et.  al.  (Ref  3)  and 
Funk  and  Wood  (Ref  5)  to  predict  this  increased  attenuation. 


Objectives 


The  following  objectives  were  established  to 


1 


investigate  the  effects  of  turbulent  mean  flow  and 
oscillatory  signals  in  pneumatic  transmission  lines  . 

1*  To  experimentally  determine  the  small  signal 
response  of  lines  with  developed  turbulent  mean 
flow. 

2.  To  modify  existing  computer  programs  to  incorporate 
new  resistance  expressions. 

3.  To  develop  models  for  phase  shift  and  attenuation. 


II.  Theor] 


Many  previous  investigations  of  fluid-transmission 
lines  have  successfully  used  a  pneumatic-electrical  analogy 
in  their  analysis  with  pressure  and  volumetric  flow  rate 
analogous  to  voltage  and  current  respectively.  Using  this 
analogy  the  following  equations  describe  the  relationship 
between  pressure  and  volumetric  flow: 


37  =  za  =  (r  +  juil)o 


dG  _ 


-  (G  +  _'uJ  C  )  P 


where  the  complex  terms  Z  and  Y  are  defined  as  the  series 
impedance  and  shunt  admittance  of  the  line. 

Krlshnaiyer  and  Lechner  (Ref  2)  present  Z  and  Y  in 
round  lines  with  blocked  or  laminar  flow  as 


L 1"  L  P  • 
?  lL  j 


+  ;Kr  +  iiy-fr/ 

“  I  A  AC  l 


uj(y-Oyp[DG. 

[dcY  +  [dgJ^ 


(T-  O  tP  [C' 

be] Y+  [c:-  ^ 


where  the  terms  DR,  DL,  DG,  and  DC  are  given  as 


DR  =  J/8  ♦  h^/4  *  3/lShv) 


(5) 


DL  =  hv/4  -  15/(54hJ 

(6) 

DG  =  ht/2  -  l/(4ht) 

(7) 

DC  =  1/4  *  hr/2  +  7//4hr J 

(8) 

with  and  hf  defined  as  velocity  and 

distribution  parameters  given  by  equations  (9) 

temperature 

and  (10) 

Vg 

h  v  =  2 1  cu  /  to  y) 

(9) 

h-f  —  <j  hy 

(10) 

U)v  is  a  characteristic  frequency  defined  by  Nichols  (Ref  1) 

as 


Ij) 


✓ 


(ID 


These  equations  are  accurate  over  a  frequency  range  of 
0.1u>v  <  u.)  <  OO  for  blocked  lines  and  for  lines  with  fully 
developed  laminar  flow. 

Transmission  lines  can  be  described  by  a  propagation 
operator  per  unit  length  P  and  a  characteristic  Impedance 
ZQ  given  by  equations  (12)  and  (13) 

r  =  (ZY)'/z  =  OC  +  j/3  (12) 

Z0^tZ/Y)Vz  (13) 

where  Or  is  defined  as  the  attenuation  and  {3  is  the  phase 


4 


angle  per  unit  length.  The  Input  Impedance  of  a  line  is 
given  as 


(ZT  +  Zo)exp(ri)  +  izr  -  Zc)expim 
s  ~  Z°(ZT  +z0)9xpiri)  -  (Zt  -  Zpjtxpirr, 


(14) 


The  ratio  of  the  recievlng  pressure  to  the  sending 
pressure  is 


Pr  _  zzT 

~  izt  +  Zc)expfri)  *  TTf  -  Zc)expiri:  (15) 

The  gain  is  the  magnitude  of  the  complex  pressure  ratio 
given  by  equation  (15) 


1  - 


and  in  decibels 


(16) 


3dt  -  20/oS|o5  (17) 

The  phase  shift  between  the  recievlng  and  sending 
points  on  the  line  is  the  angle  formed  by  the  ratio  of  the 
imaginary  to  the  real  parts  of  the  pressure  ratio 

lrr(Pr/Ps) 

Beta  -  ran'  - 2_  (18) 

Re(Pr /Pfc  ) 

Turbulent  Flow  Effects 


In  order  to  predict  the  attenuation  for  lines  with 
fully  developed  turbulent  mean  flow,  a  constant  LRC  model  is 


used  following  the  method  of  Moore  (Ref  4).  This  model  uses 
the  turbulent  resistance  of  the  line,  the  adiabatic 
lnertance  of  the  line,  and  the  isothermal  capacitance 
(adiabatic  capacitance  with  7  =  1.0  )  of  the  line  to 
determine  the  propagation  operator  and  characteristic 
impedance.  The  turbulent  resistance  is  defined  as 


(23) 


(25) 


The  turbulent  Impedance  is  a  complex  expression 
consisting  of  a  real  term,  the  turbulent  resistance,  and  an 
imaginary  term,  the  adiabatic  lnertance. 

Z  T  =  Rvi  *  JiOLe  (26) 

Using  the  adiabatic  lnertance  introduces  a  small  error  into 
the  model  but  greatly  simplifies  it  and  reduces  computation 
time.  The  actual  turbulent  lnertance  is  slightly  greater 
than  the  adiabatic  lnertance,  Moore  (Ref  4)  shows  that  the 


turbulent  lnertance  is  approximately  3.5Z  larger  than  the 


adiabatic  lnertance  at  a  Reynolds  number  of  5 ,000.  The 
turbulent  lnertance  approaches  the  adiabatic  lnertance  as 
the  Reynolds  number  increases* 

The  turbulent  admittance  is  a  complex  expression 
consisting  of  a  conductance  (real  part)  and  a  capacitance 
(imaginary  part).  The  capacitance  is  thought  to  be 
isothermal  at  low  frequencies  so  the  expression  for  the 
adiabatic  capacitance  is  used  (Eq.  (25))  with  gamma  ■  1. 
The  conductance  is  assumed  to  be  zero  for  the  low 
frequencies  at  which  this  model  is  used 

YT  =  j^Ca  (27) 

Moore  (Ref  4)  uses  the  Blaslus  relation,  equation  (28), 
as  a  simple  method  of  relating  fRe  to  the  Reynolds  number 
for  turbulent  flow. 


fRe  -  0.3164  Re3/4 


(28) 


For  laminar  flow  the  expression  relating  fRe  to 
Reynolds  number  is  Eqn.  (29). 


fRe  -  64 


(29) 


The  turbulent  resistance  then  becomes 


Rvf 


0.3164  Re^4  ^ 
2AD2 


(30 


and  the  laminar  resistance  for  this  model  becomes 


Rvl 


(31) 


The  resistance  used  in  this  model  is  a  DC  resistance, 
Franke  (Ref  6)  suggested  that  using  the  DC  resistance  might 
cause  the  constant  LRC  model  to  predict  lower  attenuations 
than  would  actually  be  encountered.  The  constant  LRC  was 
modified  to  account  for  a  possible  higher  AC  resistance. 
The  AC  resistance,  as  suggested  by  Franke  (Ref  6),  is 


Rvr^AC 


AD 


f  * 


a  f 
2  dU 


(32) 


After  carrying  out  the  differentiation  and  simplifying  th< 
expression,  the  AC  resistance  becomes 


pvr,«:  =  '•  (33) 

while  the  expression  for  the  laminar  resistance  will  not 
change . 

The  attenuation  predicted  by  the  constant  LRC  model, 
using  both  the  DC  and  AC  resistance,  is  nondimensionalized 
using  the  method  of  Brown  (Ref  3)  and  plotted  against  a 
nondlmens  lonal  freque^  Cl  in  Fig  1.  The  attenuation 
predicted  by  Nichols  equations,  as  approximated  by 
Krlshnaiyer  and  Lechner  (Ref  2),  is  also  shown  as  a 
reference.  Fig.  1  shows  that  the  constant  LRC  models  predict 
attenuations  higher  than  those  predicted  by  Nichols  below 


(Ref  7)  produced  large  deviations  In  attenuation  In  the 
vicinity  of  the  break  frequency. 

The  computer  program  used  in  this  study  employed  the 
turbulent  constant  LRC  model  at  frequencies  below  the  break 
frequency,  u>£,  ,and  the  laminar  model  (Eqns.  (3)  and  (4))  at 
frequencies  greater  than  the  break  frequency.  The  program 
uses  these  models  to  determine  the  propagation  operator  and 
characteristic  impedance  (Eqns.  (12)  and  (13))  then  uses 
Eqns.  (14)  and  (15)  to  determine  the  theoretical  pressure 
ratio.  The  program  then  uses  Eqns.  (16),  (17),  and  (18)  to 
compute  the  theoretical  gain  and  phase  shift. 

Mean  Pressure  Drop 

Schlichting  expresses  the  relation  between  pressure 
drop  and  distance  in  a  constant  area  duct  as 

A  P  =  -f  —  (34) 

D  2 

If  equation  (28)  is  used  to  determine  f  and  substituted  into 
equation  (34)  then  the  pressure  drop  can  be  expressed  as  a 
function  of  velocity  and  Reynolds  number.  Further 
manipulation  with  the  velocity  and  density  can  result  in 
equations  (35)  and  (36)  which  express  the  pressure  drop  as  a 
function  of  1,  D,  Re,  p,  and  fd  for  turbulent  and  laminar 
flows  respectively. 


Equations  (35)  and  (36)  are  only  applicable  for 
incoapresaible  flow.  To  find  the  pressure  drop  in  small 
diameter  lines  at  high  Reynolds  numbers,  Eqns.  (35)  and  (36) 
are  rewritten  as  lnfinltesmal  pressure  drops  for  a  small 
distance  dx.  The  ideal  gas  law  is  used  to  express  density 
in  terms  of  pressure  and  temperature,  then  assuming  that 
Reynolds  number  and  temperature  are  constant  in  the  line  the 
equations  can  be  Integrated  to  get  Eqns.  (37)  and  (38). 


III.  Experimental  Apparatus 


The  experimental  apparatus  consisted  of  pneumatic 
signal  generating  equipment  and  signal  analysis 
Instrumentation  shown  schematically  In  Fig.  2.  A  rotometer 
was  also  used  to  verify  that  the  mean  flow  was  consistent 
with  the  pressure  drop. 

The  wave  analyzer  connected  to  the  sending  dynamic 
pressure  transducer  was  used  to  generate  a  sinusoidal 
signal.  This  signal  was  then  amplified  and  sent  to  the 
pneumatic  signal  generating  equipment  which  amplified  it 
further  and  sent  it  to  the  pneumatic  signal  generator.  A 
frequency  counter  was  used  to  accurately  determine  the 
actual  frequency  of  the  signal  generated  by  the  wave 
analyze r . 

The  signals  in  the  line  were  measured  by  two  dynamic 
pressure  transducers.  These  signals  were  sent  to  charge 
amplifiers.  The  amplified  signal  was  then  passed  to  the 
dual  beam  oscilloscope  and  to  the  two  wave  analyzers.  The 
oscilloscope  was  used  to  monitor  the  signal  and  to  measure 
the  phase  delay.  The  wave  analyzers  were  used  to  measure 
the  RMS  output  voltage  of  the  signals. 

The  mean  flow  rate  was  established  using  the  rotometer. 
Once  the  mean  flow  rate  was  established  the  rotometer  was 
removed  and  the  flow  rate  was  monitored  using  the  static 
pressure  transducer.  The  losses  associated  with  the 


rotometer  were  small  enough  that  its  removal  caused 


negligible  changes  in  flow  rate.  The  rotometer 
calibrated  before  any  tests  were  conducted. 

Test  Configurations 


wa  8 


Six  different  line  configurations  were  used  for  the 
experiment.  A  schematic  of  the  basic  line  configurations  Is 
shown  In  Fig  3.  Each  test  line  consisted  of  plexiglas 
fixtures  at  each  end  which  contained  the  transducer  cavities 
and  a  smooth  steel  line  of  constant  diameter  between  the  two 


plexiglas 

blocks  . 

The  lines 

were 

smooth 

enough 

t  0 

accurately 

predict 

the  pressure 

drop 

in  the 

line 

using 

equation 

(37)  for 

turbulent  flow 

and 

equation 

(38) 

for 

laminar  flow. 

Classification  of  Cases 

The  cases  were  designated  using  a  two  digit  system. 
The  first  digit  indicates  the  line  being  tested  and 
corresponds  to  the  line  numbers  used  in  Table  I.  The  second 
digit  refers  to  the  Reynolds  number  at  which  the  line  is 
being  tested.  A  one  indicates  a  Reynolds  number  of  2000,  a 
two  indicates  a  Reynolds  number  of  5000,  and  a  3  indicates  a 
Reynolds  number  of  10000.  For  example  case  52  is  line  5 
tested  at  a  Reynolds  number  of  5000. 

Line  1  was  tested  in  several  different  configurations 
before  taking  data  for  all  of  the  lines  at  the  various 
Reynolds  numbers.  One  test  was  to  verify  that  the  length  of 
the  line  from  the  pneumatic  signal  generator  to  the  test 
lines  did  not  have  an  adverse  effect  on  the  signals  and  the 


Figure  3  Test  Line  Configuration 


other  to  determine  the  effect  of  the  volume  immediately 
beneath  the  dynamic  transducers.  These  two  cases  were 
tested  at  Reynolds  number  2000  on  line  1  and  were  designated 
caaea  11A  and  11B  respectively. 


IV.  Experimental  Procedure 


The  same  procedure  was  used  consistently  to  collect 
data  for  each  test.  Equation  (37)  or  (38)  was  used  to 
calculate  a  pressure  drop  for  the  line  at  the  Reynolds 
number  desired.  The  output  of  the  static  pressure 
transducer  was  then  used  to  set  up  the  equivalent  pressure 
drop.  The  rotometer  was  then  used  to  verify  that  the  flow 
rate  was  correct  for  the  case.  The  rotometer  was  then 
removed  and  the  static  transducer  was  used  to  monitor  the 
mean  pressure  for  any  variations  in  flow  rate.  Test  runs 
were  made  at  Reynolds  numbers  of  2000,  5000,  and  10000  for 
the  four  smaller  diameter  lines.  The  two  0.195  diameter 
lines  were  run  at  Reynolds  numbers  of  2000  and  5000  due  to 
the  mass  flow  constraints  of  the  pneumatic  signal  generator. 

The  barometric  pressure  and  temperature  in  the  room 
were  recorded  before  each  test  run.  The  desired  flow  rate 
was  set  using  the  ptocedure  outlined  above  then  the  wave 
analyzer  was  used  to  set  the  frequency  of  the  signal.  The 
RMS  output  of  the  sending  and  recieving  transducers  was 
recorded  from  the  wave  analyzers  and  the  phase  lag  was 
measured  on  the  oscilloscope. 


V.  Results  and  Discussion 


The  computer  program  used  to  predict  the  theoretical 
gain  and  phase  shift  was  based  on  that  used  by  Malanowskl 
(Ref  9).  It  was  extensively  modified  to  incorporate  the 
constant  LRC  turbulent  models.  In  order  to  verify  that  the 
original  laminar  section  still  functioned,  the  program  was 
used  to  predict  the  gain  and  phase  shift  in  a  blocked  line. 
Then,  the  program  was  used  to  predict  the  gain  and  phase 
shift  for  each  of  the  six  lines  with  laminar  mean  flow. 


Table  II  gives 

a  summary  of 

test 

conditions  and 

some 

experimental  and 

theoretical 

data 

for 

the 

cases 

with 

Reynolds  numbers 

of  2000. 

Figures 

4-18 

show 

the 

experimental  and 

theoretical 

gains 

and 

phase 

shifts 

for 

these  cases. 

There  are  no 

phase  diagrams  for 

line 

6  at 

any  of 

the 

Reynolds  numbers  tested;  this  is  due  to  the  fact  that  the 
noise  in  line  6  made  accurate  collection  of  phase  lag 
impossible.  This  is  also  true  of  several  of  the  other  lines 
at  higher  Reynolds  numbers. 

Figures  4-18  verify  that  Nichols  equations  are  very 
good  at  predicting  both  the  gain  and  phase  shift  for  lines 
with  laminar  mean  flow.  Only  cases  31  and  61  show  major 
discrepancies,  occuring  mainly  at  lower  frequencies.  This 
may  be  partially  due  to  noise  problems,  described  above, 
which  were  associated  with  the  two  smallest  lines.  All  of 
the  laminar  curves  however,  showed  a  1-2  db  error  at 
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frequencies  below  approximately  which  approaches  the  low 

frequency  limit  for  which  Krlshnaiyer  and  Lechners  equations 


are  applicable. 

Turbulent  Flow 

Figures  19-34  show  the  gain  and  phase  shift  for 
Reynolds  numbers  from  5,000  to  10,000.  Table  III  gives  the 
value  of  the  nondimensional  frequency  (used  in  Fig.  1) 
which  corresponds  to  the  highest  tested  frequency  of  1000  hz 
for  each  case.  The  value  of  ,  which  corresponds  to  the 
intersection  of  the  laminar  and  turbulent  attenuation  curves 
is  also  found  in  Table  III  for  both  AC  and  DC  resistance. 

The  method  used  to  predict  the  gain  and  phase  shift 
employed  the  turbulent  model  up  to  the  break  frequency  . 
From  this  frequency  up  Krlshnaiyer  and  Lechners  equations 
were  used  to  predict  the  gain  and  phase  shift.  Because  of 
t hi s  the  curves  representing  the  constant  LRC  models  are 
only  plotted  up  to  their  break  frequency. 

Figures  19-26  show  the  theoretical  gain  and  phase  shift 
and  the  experimental  gain  and  phase  shift  for  four  of  the 
lines  at  a  Reynolds  number  of  5,000.  The  figures  show  that 
Nichols  theory  is  accurate  above  the  break  frequency. 
Figures  27-28  show  large  discrepancies  in  all  of  the 
theories  although  the  constant  LRC  with  AC  resistance  is 
good  at  high  frequencies.  The  constant  LRC  models  appear  to 
be  inadequate  to  predict  the  phase  shift  for  any  Reynolds 
number  tested  while  Krlshnaiyer  and  Lechners  equations 


Figure  19  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  12 
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Figure  31  Experimental  and  Theoretical  Gain  vs  Frequency  for  Case  53 
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are  accurate  for  phase  information  for  most  of  the  lines 
tested. 

Figures  29  and  31  best  demonstrate  how  Krishnaiyer  and 
Lechners  equations  predict  attenuations  which  are  too  low. 
The  constant  LRC  model  with  DC  resistance  also  predicts 
attenuations  which  are  too  low  while  using  the  AC  resistance 
results  in  much  better  predictions  in  attenuation  and  gain. 
Both  of  the  constant  LRC  models  are  shifted  on  the  frequency 
axis  in  the  figures.  The  phase  shift  predicted  is  also  very 
far  off.  This  may  be  due  to  the  assumption  made  in  the 
models  to  use  the  isothermal  capacitance.  Moore  (Ref  4) 
showed  that  the  capacitance  in  blocked  lines  is  isothermal 
at  low  frequencies  and  adiabatic  at  high  frequencies.  The 
author  knows  of  no  method  at  this  time  to  determine  the 
frequency  dependence  of  the  capacitance  for  turbulent  mean 
flows . 

Cases  33  and  63  are  similar  to  cases  32  and  62  in  that 
the  theory  does  not  predict  the  gain  very  well.  There  are 
very  large  pressure  drops  associated  with  these  cases.  The 
models  were  all  developed  for  incompressible  flow;  this  may 
be  one  of  the  reasons  the  theoretical  curves  are  off  so 
much.  In  addition,  Mach  numnber  effects  were  not  accounted 
for.  There  is  no  phase  data  available  for  these  two  cases 
because  of  the  noise  problem  discussed  earlier. 


VI.  Conclusions 


1.  The  constant  LRC  model  with  AC  resistance  can 
accurately  predict  the  attenuation  but  cannot  predict  the 
phase  shift  as  it  is  used  here.  Because  of  this  the  model 
cannot  be  used  to  predict  the  gain  because  it  shifts  the 
frequency  of  the  harmonic  peaks.  At  this  time  no  method 
exists  to  predict  the  difference  in  capacitance  due  to 
turbulent  boundary  layers.  The  changes  in  capacitance  may 
be  important  in  matching  theoretical  solutions  with  data. 

2.  The  constant  LRC  model  is  most  accurate  when  the  AC 
resistance  is  used.  Using  the  DC  resistance  causes  the 
attenuation  to  be  underestimated  and  prevents  the  model  from 
being  used  to  predict  gain. 

3.  None  of  the  models  or  equations  used  performs  well 
in  cases  with  large  pressure  drops  (compressible  flow). 

A.  The  propagation  characteristics  of  fluid 
transmission  lines  with  turbulent  flow  exhibits  laminar 
characteristics  at  high  frequencies  and  increased 


attenuations  at  low  frequencies 


VII.  Recommendations 


1.  Studies  of  the  effects  of  turbulent  boundary  layers 
on  capacitance  should  be  conducted.  If  possible  simple 
models  should  be  developed  to  predict  the  capacitance. 

2.  Studies  on  the  effects  of  compressible  flow  on  the 
propagation  characteristics  of  fluid  transmission  lines 
should  be  conducted. 

3.  Experimental  and  analytical  research  should  be 
conducted  on  the  low  frequency  characteristics  of 
transmission  lines  with  turbulent  mean  flow. 
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ram 


Read  ICAS:  Reads  the  experiment  case  number. 

Read  TF,  PG,  AMU,  RE,  GAM,  SIG: 

TF  is  the  line  temperature  in  deg  F 

PG  is  the  atmospheric  pressure 

AMU  is  the  dynamic  viscosity 

RE  is  the  gas  constant 

GAM  is  the  ratio  of  specific  heats 

SIG  is  the  square  root  of  the  Prandtl 

number . 

Read  N:  N  is  the  number  of  lines. 

Read  DI(I),  AD(I),  P(I): 

DI  is  the  line  diameter 
AD  is  the  line  length 
P  is  the  line  pressure. 


Plotter  Subroutine 
Read  CASE(I):  CASE  same  as  ICAS. 

Read  NPTS,  LSMB : 

NPTS  is  the  number  of  experimental  input 
points 

LSMB  is  the  plotter  symbol  used  to 
represent  the  experimental  data 

Read  FREQ,  PS,  PR,  PHT 

FREQ  is  the  experimental  frequency  at 
which  the  data  was  recorded. 

PS  is  the  RMS  voltage  of  the  sending 
dynamic  transducer 

PR  is  the  RMS  voltage  of  the  recieving 
dynamic  transducer 

PHT  is  the  phase  lag  between  the  sending 
and  recieving  ports  in  milliseconds. 


m 


PROGRAM  MSB(INPUT, OUTPUT, TAPE 5-INPUT, TAPE 6-OUTPUT, TAPE 7, TAPE 9) 

C 

C  TRANSMISSION  LINE  STUDY/Z-IN  METHOD  WITH  BRANCH  SHORT  OUTPUT 
DIMENSION  DBT(300),DI(100),AD(100),AR(100),CV(100),RVT(100), 
lOMT(10O),CNA(10O),GMA(100),AGM(100),FN(I00),RN(I00),GN(100), 

2ALN( 100) ,CN( 100) ,AN( 100) , BTN( 100) , AMC( 100) ,AZRN( 100) , BZRN( 100) , 
3AZIN(100),BZIN(100),DC(100),DG(100),DD(100),P(100),RHO(100), 
4ANU(100) ,CA( 100) ,QTCA( 100) ,0MG(300) ,RTP(100) ,GP( 100) , BETA( 100) , 
5AZOT(100),BZOT(100),REY(100),Q(100),BET(100),TURB(100) 

DIMENSION  DB(300),OMGX(300),OMGP(300),PHASE(300),PH(300),CASE(100) 
DATA  PI/3.1415926/ 
nATA  TPI/6. 2831853/ 

C*********************************************************************** 

c* ********************************************************************** 
q* ********************************************************* ************* 
CALL  PL0TS(0. ,0. ,9) 

1  READ(5,*)ICAS 

110  READ(5,*)TF,PG,AMU,RE, GAM, SIG 
READ(5,*)  N 
C 

C  ******************************************************************** 
C  AD  IS  THE  LINE  LENGTH 

C  DI  IS  THE  LINE  DIAMETER 

C  P  IS  THE  INLET  PRESSURE 

C  A  ONE  FOR  NLINE  CR  MLINE  INDICATES  THE  LINE  IS  OPEN  AND  ANY  OTHER 
C  NUMBER  INDICATES  A  BLOCKED  END. 

C  ******************************************************************** 

c 

DO  898  I-1,N 

READ(5,*)  DI(I),AD(I),P(I) 

898  CONTINUE 
NPG-0 
ICT-1 
IND-1 
NNN-1 
M-0 

Q* ***************************************** *********************** 

C  DETERMINE  REYNOLDS  NUMBER  FOR  LINE 
TOTL-O.O 
DO  2  1-3, N 
TOTL-TOTL  4AD ( I ) 

2  CONTINUE 

RHO(N)-(PG"fP(N)  )/  (RE*(TF+460)  ) 

REYT-(PG*P(N)+.5*(P(N)**2.))*(DI(N)**3.)/32. /AMU/AMU 
REYT-REYT/TOTL/RE/ (TF+460. ) 

IF  (REYT.LT.2300.)  GOTO  3 
REYT-(REYT*32 . / . 1582 )**(4 . /7 . ) 

3  QM-REYT*P I *D I ( N ) *AMU/ 4 . 

DO  23  1-1 ,N 
PBR-P(I)4PG 
TBR-TF4460. 

RHO( I )-PBR/ (RE*TBR) 

ANU(I)-AMU/RHO(I) 


CA( I )-SQRT(PBR*GAM/RHO( I ) ) 

AR(I)«PI*DI(I)*DI(I)/4. 

CV(I)-(8.*PI*ANU(I))/AR(I) 

OMT(I)-CV(I)/(SIG*SIG) 

CNA(I)-(8.*PI*AMU)/(AR(I)*AR(I)) 

Q(I)-QM/RHO(I) 

23  CONTINUE 
Q(l)-0.0 
Q(2)-0.0 
GMN1-. 5*(GAM-1 . ) 

DO  26  I-1,N 

REY(I)-(4.*RH0(I)*Q(I))/(PI*DI(I)*AMU) 

PBR-P(I)4PG 
TEMP-GMNI / ( GAM*  PBR ) 

GMA(I)«TEMP*AR(I) 

AGM(I)-AR(I)/(GAM*PBR) 

QTCA(I)«0.25*CA(I) 

RVT(I)«.3164*AMU/2./AR(I)/DI(I)/DI(I) 

RVT(I)-RVT(I)*(REY(I)**.75) 

26  FN(I)-QTCA(I)/AD(I) 

NST-1 

DW-O. 

Y-l. 

40  DO  80  J-l.NST 
M-M+l 
Y-Y-HJW 
W-TPI*Y 
DO  27  I«1,N 
ARG-. 5*SQRT(W/CV(I) ) 

27  RN(I)-CNA(I)*(  .375-hARG-K  .375/(4. *ARG)  )  ) 

DO  28  I-1,N 

DC(I)-.254SQRT(W/0MT(I))+.  125*SQRT(0MT(I)/W) 

DG(I)-SQRT(W/0MT(I))-.125*SQRT(0MT(I)/W) 

DD(I)-DC(I)*DC(I)+DG(I)*DG(I) 

28  GN(I)-W*(GAM-1.)*AGM(I)*DG(I)/DD(I) 

DO  29  1-1,  N 

ARG- . 5*SqR  T ( CV ( I ) /W ) 

29  ALN(I)-RHO(I)*(l.-tARG-(ARG*(15.*CV(I)/(W*64.))))/AR(I) 
TEMP-GMNI /W 
DO  30  1-1 ,N 
TURB(I)-0 

30  CN(I)-AGM(I)*(1.+((GAM-I.)*DC(I)/DD(I))) 

TEMP— W*W 
DO  31  1-1 ,N 

TEM1-RN( I )  *GN(  I  )+TEMP*ALN(  I )  *CN(  I ) 
TEM2Hf*(RN(I)*CN(I)+GN(I)*ALN(I)) 

CALL  RTCMP(ARG1,ARG2,TEM1,TEM2) 

AN(I)-ARG1 

BTN(I)-ARG2 

TEM1— W*W*RH0(I)/(P(I)4PG) 

TEM2Hf*AR(I)*RVT(I)/  (P(I  )-fPG) 

CALL  RTCMP (ARG1 , ARG2 , TEM 1 , TEM2 ) 

TAN-ARG1 


TBTN-ARG2 

IF(AN(I) .GE.TAN)  GO  TO  31 

AN(I)-TAN 

BTN(I)-TBTN 

TURB(I)-1 

31  AHC(I)-TPI/BTN(I) 

C  CALCULATE  ZO 
DO  32  I-1,N 
TEM1-W*ALN(I) 

TEM2Hf*CN(I) 

TEM3-RN(I) 

TEM4-GN(I) 

CALL  CMPDV ( ARG1 , ARG2 , TEM3 , TEM 1 , TEM4 , TEM2 ) 

CALL  RTCMP(AZRNI,BZRNI,ARG1,ARG2) 

AZRN(I)-AZRNI 

BZRN(I)-BZRNI 

IF(TURB(I).EQ.O)  GO  TO  32 

TEMl-RHO(I)*(P(I)-H>G)/AR(I)/AR(I) 

TEM2— RVT(I)*(P(I)+PG)/W/AR(I) 

CALL  RTCMP(ARG1 , ARG2 , TEM1 , TEM2 ) 

AZRN(I)-ARG1 
BZRN(I)-ARG2 
32  CONTINUE 
C  CALCULATE  Z  IN  1 
1-1 

TEMP-AN(I)*AD(I) 

IF(TEMP.GT.88. )  GO  TO  80 
ARGl-COSH(TEMP) 

ARG2-SINH(TEMP) 

TEM5-BTN(I)*AD(I) 

TEMP-COS (TEM5) 

TEM1-ARG1*TEMP 
TEM3-ARG2  *TEMP 
TEMP-SIN(TEM5) 

TEM2-ARG2*TEMP 
TEM4-ARG 1 *TEMP 

CALL  CMPDV (ARG 1 , ARG2 , TEM 1 , TEM  2 , TEM  3 , TEM4 ) 

TEMl-AZRN(I) 

TEM2-BZRN(I) 

CALL  CMPMP(TEM3,TEM4,TEM1,TEM2,ARG1,ARG2) 

AZIN(I)-TEM3 
BZIN(I)-TEM4 
C  CALCULATE  Z  IN  3 
1-3 

TEMl-AZRN(I) 

TEM2-BZRN(I) 

TEM3-0. 

TEM4-0. 

TEM5-AN(I) 

TEM6-AD(I) 

TEM7-BTN(I) 

CALL  CALZ IN(AARG , BARG , TEM 1 , TEM 2 , TEM3 , TEM 4 ,TEM5, TEM6 , TEM 7 ) 
AZIN(I)-AARG 


62 


r%v 


BZIN(I)-BARG 
C  CALCULATE  Z  IN  4,5 
DO  34  1-4,5 
TEMl-AZRN(I) 

TEM2-BZRN(I) 

TEM3-AZIN(I-1) 

TEM4-BZIN(I-1) 

TEM5-AN(I) 

TEM6-AD(I) 

TEM7-BTN(I) 

CALL  CALZ IN( AARG , BARG , TEM 1 , TEM2 , TEM3 , TEM4 , TEM5 , TEM6 , TEM7 ) 
AZIN(I)-AARG 
34  BZIN(I)-BARG 
C  CALCULATE  Z  IN  2 
1-2 

TEMI-AZRN(I) 

TEM2-BZRN(I) 

TEM3-AZIN(I-1 ) 

TEM4-BZIN( 1-1 ) 

TEM5-AN(I) 

TEM6-AD(I) 

TEM7-BTN(I) 

CALL  CALZ IN (AARG , BARG , TEM1 , TEM2 , TEM3 , TEM4 , TEM5 , TEM6 , TEM7 ) 

AZIN(I)-AARG 

BZIN(I)-BARG 

C  CALCULATE  RECEIVING  Z  FOR  LINE  6 
1-6 

TEMl-AZIN(I-l) 

TEM2-BZIN(I-1) 

TEM3-AZIN(I-4) 

TEM4-BZIN(I-4) 

CALL  ZEBRA (AZOTI ,  BZOTI ,  TEMl  ,TEM2 ,  TEM 3 ,  TEM4  ) 

AZOT(I)-AZOTI 

BZOT(I)-BZOTI 

C  CALCULATE  Z  IN  6  INCLUDING  R -TRANSDUCER 
TEMl-AZRN(I) 

TEM2-BZRN(I) 

TEM3-AN(I) 

TEH4-AD(I) 

TEM5-BTN(I) 

CALL  CALZ IN(AZ INI, BZ INI, TEMl, TEM2, AZOTI, BZOTI, TEM3,TEM4,TEM5) 

AZIN(I)-AZINI 

BZIN(I)-BZINI 

C  CALCULATE  Z  IN  7,  8,  9,  AND  10  INCLUDING  R-TRANSDUCERS 
DO  39  1-7,10 
TEML-AZRN(I) 

TEM2-BZRN(I) 

TEM3-AZ IN( 1-1 ) 

TEM4-BZIN(I-1) 

TEM5-AN(I) 

TEM6-AD(I) 

TEM7-BTN(I) 

CALL  CALZ IN(AZ INI , BZ INI , TEMl , TEM2 , TEM 3 , TEM4 , TEM 5 , TEM 6 , TEM7) 
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AZIN(I)-AZINI 
39  BZIN(I)-BZINI 

C  CALCULATE  P5/P6  INCLUDING  R-TRANSDUCER 
1-6 

TEMP-BTN(I)*AD(I) 

CSBIl-COS(TEMP) 

SNBIl-SIN(TEMP) 

TEMP-AN(I)*AD(I) 

ARGl-COSH(TEMP) 

ARG2-SINH(TEMP) 

TEM1-ARG1*CSBI1 

TEM2-ARG2*SNBI1 

AZOTI-AZOT(I) 

BZOTI-BZOT(I) 

AZINI-AZIN(I) 

BZINI-BZIN(I) 

CALL  CMPDV(TEM7,TEM8,AZOTI , BZOTI,AZINI, BZINI) 

CALL  CMPMP(TEM3 , TEM4 , TEM7 , TEM8 , TEM1 , TEM2 ) 

TEM5-ARG2  *C  SBI 1 
TEM6-ARG1*SNBI1 
TEM9-AZRN(I) 

TEMIO— BZRN(I) 

CALL  CMPDV (TEM7 , TEM8 , AZOTI , BZOTI , TEM9, TEMIO) 

CALL  CMPMP(TEM1 , TEM2 , TEM7 , TEM8 , TEM5 , TEM6) 

TEM I -TEM 3-TEM 1 
TEM2-TEM4 -TEM2 
TEMP-TEM1*TEM1 +TEM2*TEM2 
BETAC I )-ATAN2 (TEM2 , TEM1 ) 

BET(I )-( 180 . /PI )*BETA( I ) 

RTP(I-l)-SqRT(TEMP) 

GP(I-1)-20.*AL0G10(RTP(I-1)) 

C  CALCULATE  P6/P7,  P7/P8,  P8/P9,  AND  P9/P10  INCLUDING  R-TRANSDUCER 
DO  43  1-7,10 
TEMP-BTN(I)*AD(I) 

CSBIl-COS(TEMP) 

SNBIl-SIN(TEMP) 

TEMP-AN ( I )  *AD ( I ) 

ARGl-COSH(TEMP) 

ARG2-SINH(TEMP) 

TEM1-ARG1*CSBI1 

TEM2-ARG2*SNBII 

AZINIl-AZIN(I-l) 

BZINIl-BZIN(I-l) 

AZINI-AZIN(I) 

BZINI— BZIN(I) 

CALL  CMPDV<TEM7,TEM8,AZINI1,BZINI1,AZINI, BZINI) 

CALL  CMPMP(TEM3,TEM4,TEM7,TEM8,TEM1,TEM2) 

TEM5-ARG2*CSBI1 

TEM6-ARGI*SNBI1 

TEM9^AZRN(I) 

TEMIO-BZRN(I) 

CALL  CMPDV(TEM7,TEM8,AZINI1,BZINI1,TEM9, TEMIO) 

CALL  CMPMP(TEM1,TEM2,TEM7,TEM8,TEM5,TEM6) 


TEM1-TEM3-TEM1 
TEM2-TEM4-TEM2 
TEMP-TEM 1*TEM 1  +TEM2*TEM2 
BETA( I ) -A TAN 2 (TEM2 , TEM1 ) 

BET(I)«(180./PI)*-BETA(I) 

RTP(I-1)-SQRT(TEMP) 

43  GP( 1-1 )“2  0. *AL0G10 (RTP ( I— 1 ) ) 

RTPT-RTP(5)*RTP(6)*RTP(7)*RTP(8)*RTP(9) 

GPT-20.*AL0G10(RTPT) 

BETAR-BETA( 6 )+BETA(7 )4BETA(8 )4BETA( 9 )-ffiETA( 10) 

35  IF  (BEXAR. LE.O.)  GOTO  36 
BE TAR “BEXAR -TPI 

GOTO  35 

36  IF  (BETAR.GT.-TPI)  GOTO  37 
BETAR-BETAR+TPI 

GOTO  36 

37  BETAD“(180./PI)*BETAR 
WRITE (6 , 700) Y , BETAD  ,R  TPT , GPT 

700  FORMAT (2X.6HFREQ  -,F7.0,3X,11HBETA(DEG)  “,F7.2,5X,6HGAIN 
1 1PE 12 . 4 ,  IX , 2H0R , IX , 1PE 1 2 . 4 , 2HDB ) 

DBT(M)“GPT 
(»1G(M)“Y 
PHASE (M)-BETAD 
IF(NNN.GT.l)  GO  TO  80 
WRITE(6,810) 

810  F0RMAT(1H0,13X,6HLENGTH,9X,8HDIAMETER,7X,8HPRESSURE,9X,5HC  ADB , 10X 
1.3HFN  ,  7X,7HDENSITY,5X, 11HREYN0LDS  NO) 

802  FORMAT(1HO, IX, 5HLINE  ,12, 3X,  8(1PE12.4,2X) ) 

DO  803  I-1,N 

803  WRITE(6, 802)1, AD(I),DI(I),P(I),CA(I),FN(I),RHO(I),REY(I) 

NNN-2 

80  CONTINUE 

GO  TO  (85,611),IND 
85  Y-0. 

DW-5. 

NST-200 
IND-2 
GO  TO  40 
C 

Q  ******************************************************************** 

C  ******************************************************************** 

c  ******************************************************************** 
c  ******************************************************************** 
c  ******************************************************************** 
c  ******************************************************************** 
c  ******************************************************************** 
C  THIS  SECTION  IS  FOR  THE  CALCOMP  PLOTTER.  ONE  SET  OF 

C  EXPERIMENTAL  DATA  MUST  BE  INCLUDED  WITH  EACH  RUN. 

C 

C  NPTS  IS  THE  NUMBER  OF  EXPERIMENTAL  POINTS  TO  BE  INPUT. 

C  LSMB  IS  THE  PLOTTER  SYMBOL  TO  BE  USED  FOR  GPX  AND  PHD. 

C  SEE  A  CALCOMP  PLOTTER  USERS  MANUAL  FOR  DESCRIPTION. 

C  FREQ, PS, PR,  AND  PHD  ARE  THE  INPUT  EXPERIMENTAL  VALUES. 
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c 

(J*  ********* 

c 

611  IF( ICAS)502 , 502 , 610 
610  CALL  FACT0R(0.625) 

0MG(M+1)«0. 

0MG(M+2)-100. 

DBT(M+1)—  30. 

DBT(M+2)-5.0 

CALL  AXIS(0. ,0. , 17HFREQUENCY  (HERTZ ) ,-17, 10.0,0. ,0MG(M+1 ) , 
•40MG(M+2)> 

CALL  AX1S(0. ,0. , 15HGAIN  (DECIBELS), 15, 8. ,90. ,DBT(M+1), 
«BT(M+2)) 

CALL  LINE(0MG,DBT,M, 1,0,4) 

204  FORMAT ( 6X , 1A2 ) 

READ(7,204)CASE(1) 

READ(7,*)NPTS, LSMB 
WRITE (6, 207) 

207  F0RMAT( 1H0 , 1 OX , 1HN , 4X , 4HFR  EQ , 1 OX , 2H  PS , 1 OX , 2HPR , 1 OX , 4HGAIN , 1 OX , 

•ttHPHASE) 

J-0 

DO  69  I-l.NPTS 

612  READ(7,*)FREQ,PS,PR,PHT 
GPX-2 0 . *AL0G1 0 ( PR /PS ) 

PHD— PHT*FREQ* .  3  6 

WRITE(6, 209)1, FREQ, PS, PR, GPX, PHD 
209  F0RMAT(7X,I5,5F12.5) 

DB(I)-GPX 
0MGX(I)-FREQ 
IF(PHD.GT.O)  GOTO  69 
J-J+l 
PH(J)-PHD 
OMGP(J)-FREQ 
JMAX-J 
69  CONTINUE 

QMGX(NPTS+1  )-OMG(M+l ) 

DB(NPTS-H  )-DBT(M+l ) 

0MGX(NPTS+2  )-0MG(M+2  ) 

DB(NPTS+2  )«DBT(M+2  ) 

CALL  LINE (0MGX,DB,NPTS, 1,-1, LSMB) 

CALL  PL0T(15.0,0,-3) 

CALL  FACTOR (0.62 5) 

QMG(M+l)-0. 

OMG(M+2)-100. 

PHASE(M+1)— 360. 

PHASE  (M+2  )»4  5 . 

CALL  AXIS(0. ,0. , 17HFREQUENCY  (HERTZ) ,-17, 10. ,0. ,0MG(M+1), 
■KMG(M42)) 

CALL  AXIS(0. ,0. ,21HPHASE  ANGLE  (DEGREES) ,21,8. ,90. ,PHASE(M+1) , 
4PHASE(M+2)) 

CALL  LINE ((MG, PHASE, M, 1,0, 4) 

OMGP(  JMAX+1  )-OMG(M+l ) 

0MGP(JMAX+2  )«0MG(M+2  ) 


PH(JMAX+1)-PHASE(M+1) 

PH(  JMAX+2  )-PHASE(M+2  ) 

CALL  LINE (OMGP, PH, JMAX, 1,-1 ,LSMB) 

CALL  SYMBOL(15. ,5. , .165, 5HCASE  ,0.,5) 

CALL  SYMBOL(16.,5.,.165,CASE(1),0.,2) 

1003  CALL  PLOT( 10.0, 0,-3) 

C  ******************************************************************** 

q  ******************************************************************** 

c  ******************************************************************** 

C 

501  CALL  PLOTE 

502  STOP 
END 

C  CALCULATES  ROOT  OF  A  COMPLEX  NUMBER 
SUBROUTINE  RTCMP(X,Y,A, B) 

CALL  ANGL(TEMP,A, B) 

TEMP-. 5*TEMP 

Y-A*A4B*B 

X-SCRT(Y) 

X-SqRT(X) 

Y-X*SIN(TEMP) 

X-X*COS(TEMP) 

RETURN 

END 

C  MULTIPLIES  TWO  COMPLEX  NUMBERS 

SUBROUTINE  CMPMP(X, Y,A1 ,A2, Bl , B2) 

X-A1*B1-A2*B2 

Y-A1*B2-»A2*B1 

RETURN 

END 

C  FINDS  THE  QUOTIENT  OF  TWO  COMPLEX  NUMBERS 
SUBROUTINE  CMPDV(C1,C2,A1,A2,B1,B2) 

TEMP-B1*BHB2*B2 

C1-A1*BHA2*B2 

Cl— Cl/TEMP 

C2-B1*A2-A1*B2 

C2-C2/TEMP 

RETURN 

END 

SUBROUTINE  HSINX(ARG.X) 

A-EXP(X) 

B-EXP(-X) 

A-A-B 
ARG-. 5*A 
RETURN 
END 

SUBROUTINE  HCOSX(ARG.X) 

A-EXP(X) 

B-EXP(-X) 

A-A-ffi 
ARG-. 5*A 
RETURN 
END 


SUBROUTINE  ANGL(C,A, B) 

DATA  PI/3.1415926/ 

C-ABS(B/A) 

C-*TAN(C) 

IF  (A.GT.O.)  GO  TO  5 
IA-1 
GO  TO  7 
5  IA-0 

7  IF  (B.GT.O.)  GO  TO  10 

IB-2 

GO  TO  15 
10  IB-0 

15  IA-IA+IB+1 

GO  TO  (35, 30,25, 20), IA 
20  C-C-PI 

GO  TO  35 
25  C— C 

GO  TO  35 
30  C-PI-C 

35  RETURN 

END 

SUBROUTINE  CALZ I  N(  AZ  INI ,  BZ  INI ,  AZRN1 ,  BZRN1 ,  AZ IN2 ,  BZ IN2 ,  AN1 ,  DI 1 , 
1BTN1) 

TEMP-AN1*DI1 

CALL  HC0SX(ARG1 ,  TEMP ) 

CALL  HSINX(ARG2,TEMP) 

TEMP-BTN1*DI1 
CSBI1— COS(TEMP) 

SNBIl-SIN(TEMP) 

ZR-O. 

CALL  CMPMP(TEM1,TEM2,AZIN2,BZIN2,ARG1,ZR) 

CALL  CMPMP(TEM3 , TEM4 , AZRN1 , BZRN1 , ARG2 , ZR ) 

CALL  CHPMP(TEM5,TEM6,AZIN2,BZIN2,ARG2,ZR) 

CALL  CHPMP(TEM7,TEM8,AZRN1,BZRN1,ARG1,ZR) 

A1-TEM1+TEM3 

B1-TEM2-HTEM4 

A2-TEM5-WEM7 

B2-TEM6-WEM8 

CALL  CMPMP(TEM1,TEM2,A1,B1,CSBI1,ZR) 

CALL  CMPMP(TEM5 , TEM6 , A2 , B2 , CSBI 1 , ZR ) 

CALL  CMPMP(TEM7,TEM8,A1,B1,ZR,SNBI1) 

CALL  CNFMP(TEM3 , TEM4 , A2 , B2 , ZR , SNB 1 1 ) 

TEM3-TEM34TEM1 
TEM4-TEM4-+TEM2 
TEM7-TEM7-WEM5 
TEM8-TEM8  -WEM6 

CALL  CMPDV (TEM1 , TEM2 , TEM3 , TEM4 , TEM7 , TEM8 ) 

CALL  CMPMP(AZ INI , BZ INI , TEMl , TEM2 , AZRN1 , BZRN1 ) 

RETURN 

END 

SUBROUTINE  ZEBRA(C1,C2,A1,A2,B1,B2) 

Dl-1. 

D2-0. 
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CALL  CMPDV ( ARG 1 , ARG2 , D1 , D2 , A 1 , A2 ) 
CALL  CMPDV(ARG3,ARG4,D1,D2,B1,B2) 
ARG1-ARGHARG3 
ARG2 «ARG2 +ARG4 

CALL  CMPDV(C1 , C2 , D1 , D2 ,ARG1 , ARG2 ) 

RETURN 

END 
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